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Gravitational waves from the final stages of inspiralling binary neutron stars are expected to be 
one of the most important sources for ground-based gravitational wave detectors. The masses of 
the components are determinable from the orbital and chirp frequencies during the early part of 
the evolution, and large finite-size (tidal) effects are measurable toward the end of inspiral, but 
the gravitational wave signal is expected to be very complex at this time. Tidal effects during the 
early part of the evolution will form a very small correction, but during this phase the signal is 
relatively clean. The accumulated phase shift due to tidal corrections is characterized by a single 
quantity related to a star's tidal Love number. The Love number is sensitive, in particular, to the 
compactness parameter M/R and the star's internal structure, and its determination could provide 
an important constraint to the neutron star radius. We show that the Love number of normal 
neutron stars are much different from those of self-bound strange quark matter stars. Observations 
of the tidal signature from coalescing compact binaries could therefore provide an important, and 
possibly unique, way to distinguish self-bound strange quark stars from normal neutron stars. 

PACS numbers: 04.40.Dg, 26.60.Kp, 97.60.Jd, 95.85. Sz 



I. INTRODUCTION 

Gravitational waves from the final stages of inspiralling binary neutron stars are expected to be one of the most 
important sources for ground-based gravitational wave detectors [1| . To date, LIGO observations have only been able 
to set an upper limit to the neutron star-neutron star coalescence rate of 0.039 yr~^L^Q^ 0, where Lio is the blue 
luminosity in units of 10^" L©, which translates to about 0.075 events per year in the Milky Way. This is a thousand 
times larger than the predicted rates Nevertheless, the observed neutron star-neutron star inspiral rate from the 
universe is expected to be about 2 per day in LIGO II Q. The masses of the components will be determined to 
moderate accuracy, especially if the neutron stars arc slowly spinning, during the early part of the evolution 0, Q . 

Mass measurements from inspiralling binaries will be useful, especially in constraining the equation of state through 
limits to the neutron star maximum and minimum masses, but constraints to the radius would be much more effective 
in constraining the nuclear equation of state @. Large finite-size effects, such as mass exchange and tidal disruption, 
are measurable toward the end of inspiral , but the gravitational wave signal is expected to be very complex during 
this period. Flanagan and Hinderer @ have recently pointed out that tidal effects are also potentially measurable 
during the early part of the evolution when the waveform is relatively clean. The tidal fields induce quadrupole 
moments on the neutron stars. This response of each star to external disturbance is described by the Love number 
^2 , which is a dimcnsionlcss coefficient given by the ratio of the induced quadrupole moment Qij and the applied 
tidal field Ei-j 



2R^ 

Q.,=-k2^E,,=-\E,,, (1) 



where R is the radius of the star and G is the gravitational constant. The tidal Love number fc2, which is dimensionless, 
depends on the structure of the star and therefore on the mass and the equation of state (EOS) of dense matter. The 
quantity A is the induced quadrupole polarizability. 

Tidal effects will form a very small correction in which the accumulated phase shift can be characterized by a single 
quantity A which is a weighted average of the induced quadrupole polarizabilities for the individual stars, Ai and A2. 
Since both neutron stars have the same equation of state, the weighted average A(A^), as a function of chirp mass 
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M = mJ^^TO2^^/(mi + m2)^^^, is relatively insensitive to the mass ratio toi/to2, as is shown by Hinderer et al. [loj . 
We therefore focus on the behavior of the quadrupole polarizability A of individual stars. These are related to the 
dimensionless tidal Love number ^2 for each star by ^2 = (3/2)G'Ai?~^. The Love number k2 is sensitive to the neutron 
star equation of state, in particular to the compactness parameter M/R as shown by Damour and Nagar [ll] | and the 
overall compressibility of the equation of state. In particular, the tidal Love numbers of strange quark matter stars 
are qualitatively different from those of normal matter stars. In a fashion similar to moment of inertia measurements 
from relativistic binary pulsars [l^ . an important constraint to the neutron star radius might become possible from 
gravitational wave observations. Detection of the tidal signature from coalescing compact binaries might provide an 
important, and possibly unique, way to distinguish self-bound strange quark matter stars from normal neutron stars. 

Our paper is organized as follows. In Sec. I, a new technique for the computation of tidal Love numbers is described. 
The influence of density discontinuities and phase transitions on Love numbers is discussed in Sec. II. Results of Love 
numbers for polytropic equations of state are presented in Sec. IV. Sec. V contains results for select analytic solutions 
of Einstein's equations in spherical symmetry. Love numbers for proposed model equations of state for normal stars 
with hadronic matter and self-bound stars with strange quark matter with and without crusts are given in Sec. VI, 
wherein a comparison of results between these two distinct classes of stars are also made. In Sec VII. we discuss the 
role of a solid crust on Love numbers. Our results and conclusions are summarized in Sec. VII. Relevant parameters 
required for the computation of Love numbers for analytic solutions of Einstein's equations (discussed in Sec. V) are 
to be found in Appendix A. 



II. COMPUTATION OF TIDAL LOVE NUMBERS 



The computation of tidal Love numbers is described by Thorne and Campolattaro [l3| , Hinderer [1^ , Damour and 
Nagar [llj. We use units in which G = c = 1. In terms of the dimensionless compactness parameter /3 = M/R, the 
Love number is given by 



5 

X {2/3 (6 - 3yR + 3l3{5yB, - 8) + 2/3^ [l3 - Uyn^ + l3{3yB -2) + 2/3^(1 + yn)] ) 
+ 3(1 - 2/3)2 [2 _ + 2l3{yn - 1)} log(l - 2/3)]^^ . 

Here, yn = [rH' {r) / H{r)]r=R, where the function H{r) is the solution of the differential equation 



(2) 



H"{r)+H'{r) 



+ H{r)Q{r)^Q, 



.A(r) 



where the primes denote derivatives with respect to r, and 

\ ci{r) J 

The metric functions A(r) and v[r) for the spherical star are 



eA(.) ^ 



1 - 



2m(r) 



n -1 



(3) 



(4) 



(5) 



and Cg(r) = dp /dp is the squared sound speed. Care has to be taken in the event of a first order phase transition 
or a surface density discontinuity in the evaluation of Eq. ^ because the speed of sound vanishes. We address this 
situation in the next section. 

We note that the calculation of the tidal Love number is simplified by casting Eq. ([3]) as a first-order differential 
equation for y{r) = rH' (r)/ H{r): 



ry'{r) + y{rf + y(r)e^('') [l + ^m-'^{p{r) - p{r))\ + r'^Q{r) = , 



(6) 



so that it is necessary only to determine y^ = y{R)] the value of H{R) is irrelevant. The boundary condition for Eq. 
is y(0) = 2. 



Damour and Nagar [ll| have emphasized that the factor (1 — 2(3)^ multiplying Eq. ^ makes fc2 decrease rapidly 
with compactness /3. Additionally, we note that for small compactness parameter /?, there are severe cancellations in 
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Eq. ([3]), and it is useful to expand it in a Taylor series for /3 < 0.1: 



3+y 



VR I yl-6yR-6 g I yl+3iyl-SyR + 12 n2 . yfi+62yj;+84y|+48j;fl+36 q4 /y\ 
yR {yR+W~ 7(yR+3)^ 1^ 7(i/n+3)* ^ 



5 5y%+490y% + 1272y% + 18&iyj, + 1476yR+64:8 p5 
294 (yH+3)= '■'^ 



Note that in the Newtonian limit, /3 — > 0, we have p « p, pr << 1, and one finds 



ry\r) + y(r)2 + y(r) - 6 + = , 

ci{r) 



k.iyn) ^\(^-^\. (8) 



2 V3 + 2/J? 

Equation ^ for y must be integrated with the relativistic stellar structure, or TOV, equations: [l^ [l^ 
dp{r) [m(r) + 47rr'^p(r)] [p(r) +p(r)] dmir) 



dr rir — 2m[r)) ' c?r 

We find it convenient to employ a thermodynamic variable defined by 



47r/9(r)r^ . (9) 



dh{r) = , . , (10) 

p(r) +p(r) 

as the independent variable in place of r. A stellar model can be computed specifying the value of /i(0) at the star's 
center and integrating equations for dr/dh and dm/dh. However, since these equations are divergent at the origin 
and at the stellar surface, we employed the radial variable z = instead. One therefore has 

dz 2 ^{y/^ ~ 2m) 



dh m + Airpz^/^ ' 

dm 2 ^ ^-^ 
dh dh ' 

^^'{h),/7{h)/2 = + ye^^'^^l + 47rz{h){p{h) ~ p{h))) + z{h)Q{h) , (11) 
where Q is determined by Eq. (|4]). The behavior of y near the star's center is given by 

yih) - 2 - f^- + ^^;+^^- + ^-)/^- (/.. -h) + iih. - hf) . (12) 

Also note that y^ = y(/i = 0). 

In some cases, such as with polytropic equations of state, we found it was better to use \nh as the independent 
variable. In addition, some care has to be taken in the event that dp/dh diverges at the stellar surface, which is the 
case for polytropes if the polytropic index n < 1 . 



III. THE ROLE OF DENSITY DISCONTINUITIES AND PHASE TRANSITIONS 

As Eq. (IH) for y contains the squared adiabatic speed of sound ~ dp /dp, the solution will be altered in the case 
of phase transitions within the star, for example, between the crust and the core, or in the case of a finite surface 
density such as appears in models of strange quark stars or for a uniform density stellar model. However, in the event 
that multiple charges (e.g., electric charge and baryon number) are conserved in a phase transition, the constraint of 
global charge neutrality (two Gibb's phase rules) results in a continuous pressure versus energy density curve even if 
the phase transition is of first order. The situation of a density discontinuity was elaborated in by Damour and Nagar 
[Tl| , who showed that a large discontinuity in the energy density will greatly change the value of ^2 . 

Expressing the sound speed in the vicinity of a density discontinuity as 



dp 1 dp 
dp dp 



+ Apd{p-pd), (13) 
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where pd is the pressure at the discontinuity and Ap = p(j>d + 0) — p{pd — 0) is the energy density jump across the 
discontinuity. While solving Eqs. (jlip . this discontinuity can be taken into account by properly matching solutions 
at the point of discontinuity = r{hd): 



p{rd + f^) - Pird - f^) s „Ap 

= yyd - fcj ~ 

where e — > and p ~ m(rd) / {'iirr^/i) is the average energy density of the inner (r < 7'^) core 



yird + e) = y{rd - e) — 3- = y{hd - e) - 3—, (14) 



IV. POLYTROPIC EQUATIONS OF STATE 




FIG. 1: Contours of the dimensionless tidal Love number k2 as a function of compactness /3 — M/R and polytropic index n 
(labelled along curves) for polytropes. Contours are not shown for configurations that are hydrostatically unstable (i.e., those 
with central densities larger than that of the maximum mass). 

It is useful to evaluate tidal Love numbers for polytropic equations of state p = isTp^+^Z". Love numbers in the 
Newtonian limit for polytropes have been calculated by Brooker and OUc [T^] and Kokkotas and Schaefer In 
the Newtonian limit, it is easily observed that the values for y and ^2 are independent of the polytropic constant 
K = which scales out of Eq. However, the quadrupole polarizability A = (2/3)^2^?^, and therefore the 

gravitational wave signature, does depend on K. There exist analytic solutions for the Newtonian case for polytropes 
of indices n = and 1. In the case n = 0, an incompressible fluid, cj. = 00 and the solution inside the star which 
satisfies the boundary condition at the center is simply y(r) = 2. However, the discontinuity in the sound speed at the 
stellar surface must be taken into account. According to Eq. (|T4| . yji receives a boundary contribution AiiR^p/M = 3, 
where p is the constant energy density inside the star. Therefore, for an incompressible fluid, y^ = y{r_) — 3 = — 1 
and k2 = 3/4. 

In the case n = 1, one finds 

J3/2(Wfi) „ 7r^~9 , 15 -n' 

In the above, Ji{x) is the standard Bessel function. 

Damour and Nagar [ll| , Hinderer , Binnington and Poisson [l^ have examined relativistic polytropic equations 
of state in the case of finite compactnesss. We have repeated these calculations. For each n, the polytropic constant K 
was determined from the fiducial pressure po = 1.322 x 10"^ km~^ and po = 1.249 x 10^""^ km~^ using K = PoPo^~^^" ■ 
These values are equivalent to the pressure po ^ 1 MeV fm""^ and mass-energy density pq = 94.38 MeV fm"'' (or 
a baryon density rig = 0.1 fm^'^ for the case n — 1). These values were chosen to produce reasonable neutron star 
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FIG. 2: The dimensionless tidal Love number k2 as a function of compactness /3 = M/R and polytropic index n for polytropes. 
The polytropic index n = 0.001 for the top-most curve and in multiples of 0.1 for each succeeding curve. The thickest curve 
shows results for n = 1. 
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FIG. 3: The quantity A = (2G/3)fc2i?^ in units of km^ as a function of compactness /? — M/R for polytropes of index n. 
Contours are not shown for configurations that are hydrostatically unstable. 

radii for solar mass neutron stars. For soft EOS's, n > 1, the stellar radius decreases with increasing mass up to 
the maximum mass and the maximum mass stars are relatively lighter than for stiff EOS's, n < 1. For n < 1, the 
stellar radius generally increases with increasing mass until the maximum mass is approached. The case n = 1 is 
intermediate and has a finite radius even for a star with vanishing mass. 

The results of integrating Eq. ([6]) for these polytropic EOS's are summarized in Figs. [T] and [2] which show k2 as a 
function of /3 and n. Generally, fc2 decreases with increasing n and /3. The gravitational response is proportional to 
A — (20/3)^2^?^ and this is shown for relativistic polytropes in Figs. [3]and|4l This quantity decreases rapidly with 
increasing n, and for n > 0.5, it also decreases rapidly with the compactness parameter /3. 

We have found that the results for k2 do not significantly depend on the value K in the relativistic case by altering 




0.00 0.05 0.10 0.15 0.20 0.25 0.30 

/5 



FIG. 4: The quantity A — {2G/3)k2R' , in units of km^, as a function of compactness /3 = M/R for polytropes ranging from 
n = 0.001 (top-most curve) to 3.0 (left-most curve) in increments of 0.1. Results for the polytrope n = 1 are shown as a thick 
curve. 

our fiducial values of po or po within reasonable ranges resulting in configurations of similar dimensions to neutron 
stars. Our results are the same as those of Damour and Nagar [ll|, Hinderer [l^, Binnington and Poisson jl^ to 
within numerical accuracy. 

V. LOVE NUMBERS FOR ANALYTIC SOLUTIONS OF EINSTEIN'S EQUATIONS 



0.8 




FIG. 5: The dimensionless tidal Love number k2 as a function of compactness /? — M/R for analytic solutions (see Appendix 
A) of Einstein's equations in spherical symmetry. 

It is also useful to compute the tidal response for some of the known analytic solutions of Einstein's equations 
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in spherical symmetry. All analytical solutions are scale-free; they contain essentially two parameters, the central 
energy density Pc and compactness parameter /? = GM/Rc^. Among the useful analytic solutions we will study are 
(i) the uniform fluid sphere, (ii) the Tolman VII solution [l^, (iii) Buchdahl's solution [lO, [HI, and (iv) and (v), two 
generalizations of the Tolman IV solution (23 - l2^ . The Tolman VII and Buchdahl's solutions have vanishing surface 
energy densities and are useful approximations to realistic neutron star models. The incompressible fluid and the 
generalizations of the Tolman IV solution have finite surface densities, and the latter are reasonable approximations 
of strange quark matter stars. 

It is useful to recast Eq. pT|) in the form 



— = -2 



dw 

dy_ 

dh 



w{y/w — 2xj3) 
x/3 + a{p/ pc)w^/^ 
+ ye^ - Qe^ 
2w 



a 



dx 
dh 

A 



dw 
'dh 

P 



a p 

p 



y 



p p + p 



PcC% 



„2A 



.-A 



aw- 



,(16) 



where a = AnpcB? 



i/AI,j3 = M/R and w = jR?. Therefore, we need the quantities pjpc^ p/pc^ c^, a and 



e for each analytic equation of state. In addition, for the Tolman IV solutions, which have a finite surface density, 
the boundary contribution to ya is required. This quantity, in the present notation, is —{a/ I3){ps/ pc). The quantity 
Ps/ Pc together with the above quantities are provided in Appendix A. 

As shown in Fig. [SI the two analytic solutions that most closely resemble normal neutron stars, the Buchdahl and 
Tolman VII solutions, predict values of k2 that arc similar and which closely track the results for the n = 1 polytropc 
(of course, for /3 = 0, Buchdahl's solution and the n = 1 polytrope are identical). In contrast, the Incompressible 
and Tolman IV solutions represent a significantly different family, and, as we will see, are good approximations to 
strange quark matter stars. It is clear that the two families of analytic solutions have different behaviors, and this 
foreshadows the results for the equation of state models we discuss below. Because of the scale-free character of these 
solutions, we have not shown results for A, which will scale with the assumed pc (or, equivalently, M or R.) 



VI. LOVE NUMBERS FOR MODEL EQUATIONS OF STATE 



A. Hadronic Equations of State 




FIG. 6: Mass-radius diagram for the hadronic equation of states used in this paper. Filled (open) circles indicate configurations 
with M = 1.4 Mq (1.0 Mq). The EOS notation follows Lattimer and Prakash @| and Table I. 
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EQUATIONS OF STATE 



Symbol 


Reference 


Approach 


Comp. 


FP 


Friedman & Pandharipande 


Variational 


up 


WFF(l-3) 


Wiringa, Fiks & Fabrocine 


Variational 


up 


AP(l-4) 


Akmal & Pandharipande 


Variational 


np 


MS(0-3) 


Miiller & Serot 


Field Theoretical 


np 


MPA(l-2) 


Muther, Prakash & Ainsworth 


Dirac-Brueckner HF 


np 


ENG 


Engvik et al. 


Dirac-Brueckner HF 


np 


PAL(l-6) 


Prakash, Ainsworth & Lattimer 


Schematic Potential 


np 


GM(l-3) 


Glendenning & Moszkowski 


Field Theoretical 


npH 


GS(l-2) 


Glendenning & Schaffner-Bielich 


Field Theoretical 


npK 


PCL(l-2) 


Prakash, Cooke & Lattimer [25] 


Field Theoretical 


npHQ 


SLY4 


Douchin & Haensel [26] 


Field Theoretical 


npe 




PrakaciVi (^nnkp at T .ftttiTYipr f2Pil 


Oiiark Matter 




STE 


Steiner, Fig. [11] 


Quark Matter 


Q {u, d, s) 


PAG 


Page, Fig. m 


Quark Matter 


Q {u, d, s) 


ALF 


Alford, Fig. [TT] 


Quark Matter 


Q {u, d, s) 


HS 


Haensel, Salgado & Bonazzola [27] 


Crust 


Z,e,n 


BPS 


Baym, Pethick & Sutherland [28] 


Crust 


Z,e,n 



TABLE I: Approach refers to the underlying theoretical technique. Composition (Comp.) refers to strongly interacting 
components (n=neutron, p=proton, Z=nucleus, H=hyperon, K=kaon, Q=quark); all models include leptonic contributions. 
This table is slightly expanded from the version found in [2^ which contains references not noted here. 

The hadronic EOS's were taken from a compilation by Lattimer and Prakash @ that describes their origins. There 
are three generic families of equations of state: (i) normal nucleonic equations of state, (ii) equations of state with 
considerable softening above the nuclear saturation density, due to Bose condensation, hyperons or a mixed quark- 
hadronic phase, and iii) strange quark matter stars. We have used a selection in an attempt to span the extreme 
range of models of each type. The mass-radius curves for hadronic EOS's are shown in Fig. [6l 
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FIG. 7: The dimensionless tidal Love number k2 as a function of compactness /3 = M/R for hadronic EOSs. Filled (open) 
circles indicate configurations with M = 1.4 Mq (1.0 M0). The EOS notation follows Lattimer and Prakash and Table I. 

Love numbers as a function of compactness are shown in Fig. [7] for hadronic models. There is a relatively narrow 
spread of values of fc2 for a given compactness, and for each EOS, the value of ^2 appears to be a maximum for masses 
near 1 Mq. In contrast to the analytic Tolman VII and Buchdahl solutions, for which k2{f3 — ?• 0) ~ 0.3, fc2 tends 
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FIG. 8: The Love number k2 as a function of radius R. Filled (open) circles indicate configurations with M = 1.4 Mq (1.0 
Mq). The EOS notation follows Lattimer and Prakash 0] and Table I. 




FIG. 9: The quantity A = (2/3)fc2-R^ for hadronic equations of state. Filled (open) circles indicate configurations with M = 1.4 

Mq (1.0 Mq). 

to zero for small /3 for realistic equations of state. The fact that hadronic equations of state have a small range of 
variations as a function of compactness is reminiscent of the situation for the moment of inertia [l^ . 

It is useful to examine ^2 as a function of neutron star radius, as shown in Fig. |S1 Although the range of values 
observed for k2 are common to all models, it is now clear that the quadrupole response will vary more widely, due 
to it being proportional to R^. In Figs. and [TUl the quadrupole response is shown. The maxima in A occurs near 
1 Mq, as it did for ^2, and their is a pronounced trend for A to increase with R. Assuming the true neutron star 
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R (km) 



FIG. 10: The quantity A = {2/3)k2R'^ for hadronic equatfons of state. Filled (open) circles indicate configurations with M — 1.4 
Mq (1.0 Mo). 
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FIG. 11: Pressure versus energy density for strange quark matter equations of state with and without crust. Equation of state 
STE is taken from Steiner [301 , PAG from Page [3j| and ALF from Alford [s^l (see Table I). Density discontinuities are as 
indicated. 

equation of state is hadronic, it therefore appears that a measurement of A translates into an estimate of R relatively 
independently of the details of the equation of state. In fact, compared to the moment of inertia which scales as R^, 
the potential for a radius constraint is enhanced due the R^ behavior of A. 
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FIG. 12: Mass-radius curves for strange quark matter equations of state. The insert shows results on a logarithmic scale to 
highlight the effects of a hadronic crust. 



B. Self-bound strange quark matter stars 

Wc turn now to examine results of Love numbers for self-bound strange quark matter stars. It is uncertain whether 
or not strange quark matter stars will have signifieant crusts or not, so we examine models of both kinds. Models 
without crusts are characterized by quark matter extending up to a bare surface with a finite baryon density of 2 
to 3 times nuclear matter equilibrium density. Crusts of normal matter on top of such stars might be supported by 
stron g el ectric fields at the surface. Fig. [TT] shows three examples for both cases (STE from Steiner [s^l, PAG from 



Page [31| and ALF from Afford [32| ) . The crust and the core regions arc apparent from the large discontinuity in the 
energy density. The existence of a crust results in large radii for small stellar masses (of order 0.01 Mq), but do not 
dramatically affect the radii of stars with masses larger than 0.1 Mq (sec Fig. [T2|) . It therefore appears unlikely that 
the existence of a crust has a pronounced effect on the Love number or quadrupole properties of the star. 

In Fig. [131 the dimensionless Love number k2 is shown as a function of compactness. As was the case for hadronic 
stars, there is a clustering of curves relatively independent of the EOS for stars without crusts. The curves follow 
the analytic results for the incompressible fluid and for the Tolman IV solutions, and differ from hadronic cases by 
having a large, finite value of k2 for small /?. However, in the case of an added crust, k2 is reduced at small values of 
M/R, but this only occurs for ultra- low mass stars. For masses in excess of 1 M0, the Love number approaches the 
corresonding values for hadronic stars, and the effect of the crust is negligible. 

The quadrupole response A — 2fc2i?^/3 is shown in Fig. [14] as a function of radius. The strong dependence on 
radius follows the trend noted for hadronic stars. The effect of the crust is unimportant. 



C. Comparison of normal and self-bound stars 



In order to elaborate the distinction between strange quark matter and hadronic models, we show the quadrupole 
response A = 2fc2i?^/3 in Fig. [15] for a representative sample of models of each type. The strong dependence of 
A on i? is common to all models. Where the radii of models overlap, however, it appears that the strange quark 
matter configurations have values of A about 50% larger. This difference is probably too small to be observable, and 
it appears doubtful that any quark matter configurations will have a strong enought tidal signature to be observed. 
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FIG. 14: Quadrupole polarizabilities A for the strange quark matter stars. Filled (open) circles indicate configurations with 
M = 1.4 Mq (1.0 Mq). 
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FIG. 15: Comparison of quadrupole polarizabilities A for normal and strange quark matter stars. Filled (open) circles indicate 
configurations with Af = 1.4 Mq (1.0 Mq). 



VII. DISCUSSION 



The combined tidal effects of two neutron stars in circular orbit can be found from a weighted average of the 
quadrupole responses 

A = 4 + M)^ + (llmi + M)^) , (17) 

ZD \ mi 7712 / 

where M = mi + m2 is the total mass of the binary and Ai and A2 arc the quadrupole responses of mi and TO2. Note 
that if mi = m2, then Ai = A2 = A. If m2 = 0.5mi, then A w (40/26)Ai. It is unlikely that the mass ratio would 
be smaller than this amount, as the minimum neutron star mass that can be formed in supernovae is not less than 1 
Mq and the maximum neutron star mass is of order 2 Mq. Therefore, the value of A is similar to that of the largest 
neutron star. In the case that the individual masses can be found to reasonable accuracy from the gravitational wave 
signal, the individual values of A for the two stars will be determined to an accuracy constrained by the errors in A 
and the masses. 

We have assumed in evaluating the Love numbers that the crust behaves as a liquid. However, if the stress on the 
solid crust produced by the tidal field is large enough, then the crust can be melted and our calculations become 
valid. The strength required to melt the crust can be estimated from the results of recent work on crust breaking. 
We estimate the induced quadrupole moment to be 



3, M 



Q22 = XE22 « \^E,,E^i = y -A-^, (18) 

where the tidal field strength Eij [3^ depends on the distance D between the stars and M is the total mass; we assumed 
for simplicity an equal-mass binary. Assuming a binary in circular orbit, we can calculate the orbital frequency 17 
from Kepler's third law 

Eliminating M/ using Eq. (|T8|) , and recognizing that the frequency of the emitted gravitational waves / is twice 
the orbital frequency [3J], we have 

1/4 



27r TT V A V 3 



14 




FIG. 16: The maximum quadrupole moment Q22,max of a solid crust as a function of mass for normal and strange matter stars. 
Filled (open) circles indicate configurations with M = 1.4 (1.0 Mq). 



which has an implicit mass dependence through Q22 and A. For a 1 Mq neutron star using the EOS labelled SLY, 
Horowitz [3^ estimates that the maximum value of Q22 reached at the breaking point of the crust, where the strain 
a ~ 0.1, Horowitz and Kadau (36j . is approximately Q22,max = 10^'^ g cni'^. The breaking point is therefore reached 
during the inspiral of an equal-mass binary at the moment when the frequency of detected gravitational waves becomes 



21036gcm2 s2 

where we used the value for A for a 1 Mq star as determined in Fig. |9l Note that this frequency implies a binary 
separation distance Dtr ~ 400 km from Eq. (fT9| . Therefore, when D < D^r or / > fi,r the shear from induced 
quadrupole moment is strong enough to break the crust and beyond this point a solid crust can no longer exist. This 
frequency is below the observable region from 100 to 1000 Hz for current and proposed gravitational wave detectors 
such as LIGO Q. Consequently, during the last stages of inspiral that are observed in gravitational waves, effects 
stemming from the solid crust are probabl y ir relevant and our calculations assuming a liquid phase should be valid. 

Using the expressions provided by Owen [37[, which are supported by our results, we can approximate the maximum 
quadrupole moment for a solid crust through 

2.4 X 103« {jJL.f"' {^) neutron stars, 

^22, maa; = -jj^^ gcm^ < 3.5 X 10"^^ (sEi)^ { ^^M^ ^ hybrid and mcson-coudensate stars, , (22) 

2.8 X 10^ \,,3, ^ , , (TT^f ( solid strange stars, 

4 10^^ erg/cm vlOkm/ y M J » ' 

where amax = 0.1 is the breaking strain of the crust and yu « 4 x 10'^^ erg/cm'^ is a typical shear modulus of a strange 
quark matter crust (Horowitz and Kadau [1^), which is a thousand times the typical value in the crust of a normal 
neutron star. The results are shown in Fig. 1161 For stars with masses heavier than IMq the maximum quadrupole 
moments are within an order of magnitude of the typical value of lO^'^gcm^. 

Fig. [T7| shows results for the breaking frequency fbr calculated utilizing Eq. (|20p with the appropriate values for 
Q22,max from Fig. 1161 The breaking frequency for both kinds of stars heavier than a few tenths of a solar mass is well 
below the LIGO lower boundary of 100 Hz Q. Therefore, the crust may be assumed to be melted during the time it 
is observed, and the approximation of treating the entire star as a liquid is justified. 
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FIG. 17: The frequency of gravitational waves from an inspiraling binary when tidal forces are expected to break the crust for 
normal and strange quark matter stars. Filled (open) circles indicate configurations with M = 1.4 M0 (1.0 M0). 



The quadrupolc polarizabilities of normal neutron stars and self-bound quark matter stars have been calculated 
for a wide class of proposed equations of state of dense matter for both normal and strange quark matter stars. The 
quadrupole polarizabilities A = 2R^k2/{'SG) are characterized by the dimensionless Love number ^2 and both are 
sensitive to the equation of state, in particular to the compactness parameter M/R and the overall compressibility of 
the equation of state. For normal neutron stars, fc2 and A exhibit pronounced maxima for configurations with masses 
close to a solar mass for most equations of state. The maximum value of ^2 is not very sensitive to the EOS, lying in 
the range 0.1-0.14. In each case, maximum mass configurations have significantly lower values of fc2 and A than their 
solar mass counterparts. 

Love numbers for self-bound strange quark matter stars with or without crusts arc qualitatively different than 
those of normal neutron stars. The maxima in the value of ^2 for strange quark matter stars without crusts occurs 
for masses less than 0.1 Mq, and maximum values of order 0.8 are achieved. As in the normal matter case, the 
maxima in quadrupole polarizabilities occurs for configurations near 1 Mq. In contrast, the magnitudes of quadrupole 
polarizabilities of strange quark matter stars are usually much less than those of normal stars, owing to the larger 
radii of the latter. 

Our investigations also point the need to examine the core-crust interface region of both normal and self-bound 
quark matter stars more closely. The important issue that bears close scrutiny is the precise nature (first or second 
order) of possible phase transitions. In the case that strong discontinuities exist near the core-crust interface of strange 
quark matter stars, dimensionless Love numbers are suppressed for low mass stars relative to the cases for which there 
is no crust. However, for stars of order 1 Mq or larger, the presence or absence of a crust has little influence on Love 
numbers. 

The strength of the tidal signatures from coalescing compact binaries is proportional to A, and is therefore quite 
sensitive to the radii of the stars. For stellar configuratons with radii of order 11 km or less, the tidal response might 
be too small to observe, implying that a positive detection might be sufficient to rule out the presence of a self-bound 
star, such as a strange quark matter star, in the observed system. 



VIII. SUMMARY AND CONCLUSIONS 
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Appendix A: Parameters for Analytic Solutions of Einstein's Equations 

We use the notation /3 ~ GM/Rc^, a = iirpcR^ and x = {r/Ry. 



1. Uniform Density (p — pc) 



a = 3/3, e"^ = 1 - 2Px, 

ci^^, ^ = 1 . (Al) 



p VI - 2/3 - x/l - 2f3x 2 .„ Ps 



Pc VI -2/32;-3yi~2;8' 

2. Tolman VII (p = pc[l - x]) ^ 
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a = y/3, e-^ = 1- Px{5 - 3x) 



P 



2/3 ,1a; 

— — \ TT-r tan m 1 , 

Pc 15V/3e^ ^ 3 5' 



Wi — w 



tan 



13 



w — In 



3/3 



Wi — In 



3(1-2^)' 

1 /l - 2/3' 



6 



3/3 



tan( 



tan ( 



3e 



:;^(5-6x) 



(A2) 



3. Buchdahl's Solution (p = 12^^ - 5p) ^ 



1-/3 



smz 



1-2P ' 1-/3 + M 

A_ (1-2/3)(1-/3 + m) 



TT^X , 



^2 ' 



z (1-/3-m)(1-/? + /3cos0)' 

p _ (1-2/?)(2-2/3-3m)u p_ /3(1 - 2/3) 



1 - /3 - 4m ' 



p, (2-5/3)(l-/3 + u)2 /3' 



(l-/3 + u)2(2-5^) 



(A3) 



4. Generalized Tolman IV (N=l) (H-IH 



_ 3/3 2 - 3/3 ^ _ 1 - 3/3 + 2/3a; 

"~Tr^' " (l-3/3 + /3a;)(l-/3a;)' 

p _ 1 - 3/3 (2 - 3/3)(l - 3/3) + ^(3 - 7/3)x + 2/32a;2 _ (1 - 2/3)(l - 3^9) 

" 2 - 3/3 (1 - 3/3 + 2/3a;)2 ' ~ (1 - 3/3/2) (1 - /3) ' 

p _ 1 - 3/3 1 - X 2 _ 1 - 3/3 + 2/?x 

~ 2-3/3^1-3^ + 2/3a;' ^ 5 - 15/3 + 2/3a; ' 

5. Generalized Tolman IV (N=2) (H-IH 



a = 3/3 
^=11 

p_ 

Pc 



2-2/3 
2-5/3 

5/3x 



2/3 



3(2 - 5/3) 

2/3 



' = 1-2 



2-2/3 



2 - 5/3 + 3/3x 



2/3 



2-5/3 



5/3 



p,^(^- 5/3/3) (1-5/3/2)2/3 

Pc 



2-5/3 

2-213) 3(2 - 5/? + /?a;) 



2-2(3 



2 - 5/3 + 3/3x 



(1 - /3)5/3 
2/3 

(2 - 5/3 + 5/3a;) 



2 - 5^ + iPx 
5(2 - 5/3 + /3a;)3 



(2 - 5;3 + 3;3a;)5/3 
(2 - 2/3)2/3 



+ (2 - 5/3)^ - 5/32x2 
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